Multi-scale model of gradient evolution in turbulent flows 
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A multi-scale model for the evolution of the velocity gradient tensor in fully developed turbulence is 
proposed. The model is based on a coupling between a "Restricted Euler" dynamics [P. Vieillefosse, 
Physica A, 14, 150 (1984)] which describes gradient self-stretching, and a deterministic cascade 
model which allows for energy exchange between different scales. We show that inclusion of the 
cascade process is sufficient to regularize the well-known finite time singularity of the Restricted 
Euler dynamics. At the same time, the model retains topological and geometrical features of real 
turbulent flows: these include the alignment between vorticity and the intermediate eigenvector 
of the strain-rate tensor and the typical teardrop shape of the joint probability density between 
the two invariants, R — Q, of the gradient tensor. The model also possesses skewed, non-Gaussian 
longitudinal gradient fluctuations and the correct scaling of energy dissipation as a function of 
Reynolds number. Derivative flatness coefficients are in good agreement with experimental data. 

PACS numbers: 



The spatio-temporal fluctuations of small-scales in 
three-dimensional turbulent flows are among the most 
complex phenomena known to classical physics, being 
both highly non-Gaussian and strongly long-range cor- 
related. For example, velocity gradients as well as veloc- 
ity increments between two points typically show strong 
fluctuations much larger than their standard deviation 
P, Q ■ The same is true for fluid accelerations or velocity 
increments at two different times 0, S 0|: with long- 
range correlations up to the time scales of the largest 
eddies in the flow [B| . A possible mechanism for the large 
fluctuations is the nonlinear self-stretching [|| 0| that oc- 
curs during the Lagrangian evolution of the velocity gra- 



A systematic analysis of the dynamics of velocity gra- 
dients was made by Vieillefosse [6[. He started from the 
exact equations governing the Lagrangian evolution of 
A in the incompressible Navier-Stokes (NS) equations: 
jjjAij = —A^Akj — didjp + vd 2 Aij, where p is the pres- 
sure divided by density, v is the kinematic viscosity and 
with d/dt we mean the Lagrangian derivative. Then, he 
retained only the isotropic part of the pressure Hessian, 
didjp ~ 5ijd 2 p and used the imcompressibility condition 
to express the pressure Laplacian in terms of A. Fi- 
nally he neglected the viscous contribution arriving to 
the closed "Restricted Euler" (RE) equations: 



dients, A t 



This local self-stretching must be 



coupled with the energy exchange among larger/smaller 
eddies and with velocity fluctuations at different spatial 
locations via the pressure term. The non-linear coupling 
among different scales also relates to the concept of en- 
ergy cascade, often invoked to explain the growth of non- 
Gaussianity going from large to small scales. 
Significant advances in experimental techniques now al- 
low to measure all components of A in different spatial 
locations 0,0, E3] m high Reynolds number flows. These 
experimental measurements, together with data gener- 
ated using direct numerical simulations, has uncovered 
the existence of many interesting, and possibly univer- 
sal, geometric features of A. Namely, (i) the preferred 
alignment of the vorticity vector with the eigenvector 
of the intermediate eigenvalue of the strain-rate tensor, 
S = (A + A T )/2; (ii) the axisymmetric character of local 
deformation (two positive and one negative eigenvalue of 
S); (hi) the typical teardrop shape of the joint proba- 
bility distribution, P(R,Q) where R = -Tr(A 3 V3 and 
Q = -Tr(A 2 /2) are two invariants of A @,BIll|,tl3,[Il- 



—A = -A 2 + Tr(A 2 )I/3. 

dt 



(1) 



It is remarkable that such a simple system is already suf- 
ficient to explain many of the geometrical trends found 
in the real gradient evolution, as shown by @, 0|- On 
the other hand, the self-stretching mechanism is not con- 
strained by any energy exchange/loss mechanism, lead- 
ing to a finite time singularity for any initial condition. 
This blow up prevents the use of the RE dynamics to 
make any systematic assessment of the gradient's sta- 
tionary statistics. Prior models that seek to regularize 
the RE dynamics include a stochastic model with log- 
normal statistics of the dissipation [l5| , a linear and non- 
linear damping model for the viscous term [l6| . and a 
model where the material deformation history described 
by a tetrad of points plus some stochastic terms are used 
to mimic the anisotropic pressure fluctuations [l3j]. Re- 
cently, a model of the anisotropic pressure Hessian and 
of the viscous term has been proposed by introducing a 
finite-time memory effect in the closure of the material 
deformation tensor [T3| ■ Results show that this model re- 
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FIG. 1: Log- log plot of the mean energy dissipation (di- 
vided by 2u): (Tr(S„S^)} versus k n , for Re x = 1,500 (o); 
Re\ — 600 (□); Re\ = 130 (x) . The dissipative scale kd is 
defined as the wavenumber where the energy dissipation peaks 
(arrows in the plot). The straight line correspond to the Kol- 
mogorov spectrum ~ k^j z . The data which saturates in the 
viscous range, •, are for the coarse-grained (defined later in 
the text) variables: (Tr(S CG S CGT )). In the inset we show 
the Reynolds dependency of the energy dissipation evaluated 
at kd- The straight line corresponds to the expected slope 
oc Re\. 



produces stationarity, non-Gaussianity of transverse and 
longitudinal gradients and the correct geometrical prop- 
erties of A. Remarkably, the system also predicts exper- 
imentally observed relative scaling exponents character- 
izing how non-Gaussian statistics evolves with Reynolds 
number, at least from small to moderate values of it. 
However, attempts to simulate velocity gradient dynam- 
ics at arbitrarily large Reynolds numbers with this model 
have been, so far, unsuccessful. It was concluded that 
the difficulties were associated to the assumptions that 
the velocity gradient tensor evolves mainly uncoupled 
from larger/smaller eddies and neighboring locations. In 
this paper we propose a minimal model in which deter- 
ministic couplings among scales of motion are included. 
The model couples the self-stretching local dynamics a la 
Vieillcfossc with a cascade mechanism, i.e. by coupling 
the velocity evolution on the gradient scale given by ((T|) 
with velocity fluctuations at larger scales. The most re- 
markable result is that such a coupling will be shown to 
be sufficient to regularize the Vieillefosse finite-time sin- 
gularity, without destroying the main positive features of 
(fl]) including the presence of skewed longitudinal gradi- 
ents and long tails in the probability density functions of 
gradients; the vorticity alignment and the typical tear- 
drop shape of the joint probability, P(R,Q). Also the 
intermittency level as measured by the flatness coeffi- 
cient is in good quantitative agreement with experiments 
[H, |20T ] . The model thus provides meaningful results even 
for very large Reynolds number (here we present results 



up to Re x = 1,500). 

To introduce a cascade dynamics we start from a de- 
composition of the gradient tensor into band-passed con- 
tributions A = J2 n -^-n, each A„ describes the velocity 
gradient at a typical wavenumber k n [l8[ • The set of pos- 
sible wavenumbers is chosen equispaced on a logarithmic 
scale, k n = 2 n ko, such as to optimally capture power-law 
behaviors in the inertial range. The band-passed version 
of the NS equations is of the form 



^A 

dt " 



= -£(A p A g )„ 

Pi9 



Tr(%) 



vd 2 A n , (2) 



where the trace term is added to keep the entire right- 
hand-side trace-free, and the term B„ represents pres- 
sure effects, interscale interactions and additional spa- 
tial transport introduced by the band-pass filtering act- 
ing on the advective term on the left side of the origi- 
nal equation. Neglecting the viscous term and B„, and 
keeping only the fully diagonal term (p = q = n) in 
the double sum leads to the most severe approximation. 
Namely that the gradients, A„, band-passed on differ- 
ent shells, follow the RE dynamics separately shcll-by- 
shell: dA n /dt = -A£ + l/3Tr(A£)I. Of course this 
set of uncoupled equations suffers from the same draw- 
backs of the original Vieillefosse model. The coupling 
between different scales must be introduced by keeping 
some terms besides the purely diagonal one in the first 
term in the rhs of Eq. [2j i.e. we include a new term F 
and we choose its form by imposing an energy preserving 
structure for the coupling terms, following the assump- 
tion that the pressure-Hessian does not perform work on 
the band-passed gradient Moreover, it is natural to 
suppose that the most relevant dynamical interactions 
happen within nearby scales, (locality assumption) and 
therefore we limit the range of interactions to the next- 
nearest neighbors (up to n±2, for each k n ). One possible 
choice for the coupling term on the rtth shell is: 



¥ n [A,A]=A n+2 Al_ 
(1 - 6)2 4 A„_ 2 A„„ 1 - Tr[ 



62 2 A£_ 1 A r 



(3) 



where b is a free parameter (always fixed to b = 0.5 here- 
after). Let us notice that the structure of F n L4, A] is cho- 
sen such that the total energy E = J2 n 2 Tr(A„A^ ) is 
preserved by the effects of the three nonlinear terms. The 
model form is proposed as a model for the neglect of all 
the truncated terms in the double sum and the B„ trans- 
port term. The functional form proposed for F n [^4, A] is 
motivated by the typical structure of the nonlinear terms 
used in Shell models of turbulence (3], with the impor- 
tant difference that via ([3]) we can access now the whole 
geometrical properties of the band-passed variables, ow- 
ing to the tensorial structure of the A„ variables. To 
introduce a characteristic wavenumber from which the 
gradient variables receive their maximal contribution, we 
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FIG. 3: PDF of the cosine between the vorticity and the 
eigenvectors associated with the Intermediate (+), maximum 
(*) and smallest (x) eigenvalues of the strain-rate matrix. 
Notice the good agreement between the DNS (right) data and 
the model data (left). 



FIG. 2: Isolines of P(R, Q) for DNS and model data. Right 
panel: DNS data at Re\ ~ 180. Left panel: model data at 
Re\ = 140, evaluated at the dissipative scale kd(Re). The 
continuous line in all diagrams correspond to the Vieillefosse 
zero-discriminant curve Q 3 + 27/4R 2 — 0. 



need a viscous damping. The simplest dimensionally con- 
sistent way to do this is by modeling the viscous Lapla- 
cian term in NS equations using a linear damping term on 
each band proportional to vk^ where v is the viscosity. 
Combining all these elements one obtains the following 
dynamical system: 



dA„ 

dt 



^Tr(A^)I) 



(1 - a) (P n - vk 2 n A Tl 



(4) 

where a is a parameter which weights the relative im- 
portance of the RE dynamics with respect to the energy 
exchange and energy dissipation terms. As a first step, in 
this work we restrict attention to the case a — 0.5. The 
system consist of 9N coupled nonlinear equations. In the 
previous equations the infrared and ultraviolet trunca- 
tion is imposed by keeping in the nonlinear term A n = 
for n = {-1, -2, N+ 1,7V + 2}. 

The main question we want to answer now is if a sta- 
tionary statistical state can be achieved once a forcing is 
applied at large scales. In other words, we want to un- 
derstand if the simple energy-exchange term, F„L4, A], 
is sufficient to regularize the dynamics of the RE struc- 
ture at all scales. Hopefully, this will allow us also to 
study the important questions related to the dependency 
of geometrical structure of turbulence on the Reynolds 
numbers and its correlation with inertial range quanti- 
ties. The simplest way to achieve a stationary state is to 
add a white-in-time Gaussian variable, G, at the largest 
band, n = 0, where in order to satisfy the requirements of 
homogeneity and isotropy, we must chose the covariancc 
of the forcing to be (GijG m i) = 2S im 5ji - 1/2% <5 m ; - 



l/26u8j m [17 1 



The gradient statistics will be studied by looking at the 



shell kd (where n = d), defined as where the spectrum 
of the A n has its peak, i.e. where a balance between 
the quadratic terms and the viscous terms of (|3]) is ob- 
tained. The system of equations is integrated numer- 
ically using a 4-th order Runge-Kutta scheme. Three 
Reynolds numbers are simulated by using viscosities 
equal to v = {1CT 4 ,2.7 x 10~ 6 ,4.4 x 1CT 7 } and with 
a total number of shells N = {14,18,22}, correspond- 
ing to Re\ = {130,640, 1,500} respectively Results are 
shown in Fig [T] which shows the time-averaged energy 
dissipation (without multiplication by v), measured on 
the band-pass gradient variables, (Tr(A„A^)) as func- 
tion of k n : notice that it has the expected Kolmogorov 
scaling (Tr(A„A^)) ~ (Tr(S n S£)) ~ kt /3 . On dimen- 
sional grounds we expect that the wavenumber where 
the matching is achieved scales with Reynolds number 
as kd ~ Re\- This is indeed verified in the inset of Fig 
((T|) . Geometrical properties of small-scale turbulent fluc- 
tuations can be monitored by studying the joint statis- 
tics of the two invariants, Q,R at different Reynolds 
number and/or at different scales. Here, by measuring 
the simultaneous distribution of Q n = Tr(A^) and of 
R n = Tr(A^) for different Reynolds numbers we can 
do both. In Fig @ we show the isolines of the joint 
probability distribution P{R n , Qn) measured on the dis- 
sipative scale n = d. We notice the presence of high 
probable fluctuations along the right Vieillefosse tail, also 
seen in real turbulence, which is shown in the same figure 
panel (b) . The data are from direct numerical simulation 
(DNS), done at a similar Reynolds number. As can be 
seen, there is very good agreement of the model with the 
Navier-Stokes data, except for a small depletion of events 
in the third quadrant (i? < 0, Q > 0). In this quadrant, 
real data are strongly affected by vortex stretching, an 
effect which is evidently under-represented by the model 
evolution. The P(R, Q) distribution has very little sen- 
sitivity to Reynolds number (at least in the range inves- 
tigated here). The P(R,Q) distribution becomes more 
symmetric for band-pass variables at larger scales (not 
shown), and this effect is also observed in Navier-Stokes 
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turbulence [l2|, EH. Finally, in Fig. ([3]) we present quan- 
tification of the alignment between the vorticity at the 
dissipative scale, (wi)d = eijl(Aji)d: and the three eigen- 
vectors of the strain rate tensor also at that scale. It 
is well-known from DNS and experiments that vorticity 
tends to preferentially align with the intermediate eigen- 
vector d, 0, [lj[ : Figure (|3]) shows that our model is able 
to capture also this feature. Finally, we document the 
results in terms of intermittency. First, we show in the 
inset of Fig ([4]) the pdf of both longitudinal and trans- 
verse gradients. Correctly, the model possesses longitudi- 
nal skewed distribution while the transverse gradient pdf 
is fully symmetric. The gradient fluctuations are highly 
non-Gaussian, as a results of the growth of intermittency 
going from the large scale down to the dissipative scale. 
The growth of intermittency at decreasing length scales 
is often characterized by measuring the flatness coeffi- 
cient of velocity gradients, F 4 = {A\) / '{A d ) 2 and plot- 
ting it as function of Reynolds number Re\, where with 
Ad we mean any of the longitudinal components of A. 
As a technical difficulty we point out that the exponen- 
tial fall-off of the gradient spectrum in the shell-matrix 
model (j4]) prevents us from a direct measurement of 
on the dissipative shell (the flatness tends to grow ex- 
ponentially in the dissipative range, making the estimate 
very delicate) . For that reasons we decide to measure the 
flatness using a coarse-grained variable. Specifically, the 
coarse-grained velocity gradient at scale k m is defined as 
A^f* = Yln=o A n - For large m, we obtain a variable 
which includes fluctuations on all scales, and which be- 
comes independent on m (saturates) beyond the viscous 
range. The two variables A n and A^ G have obviously the 
same scaling in the inertial range; in the dissipative range 
the band-passed decays exponentially while the coarse- 
grained saturates, see Fig {]]). This saturation allows us 
to have robust measurements of intermittency in the dis- 
sipative range. In the main body of Fig. ((4|) we show 
the behaviour of Flatness of the coarse-grained velocity 
gradient function of Reynolds number, superposed 
with experimental data. Within the range of Reynolds 
numbers considered, the agreement is very satisfactory. 

In conclusion, we have introduced a "shell" version of 
the RE dynamics which is free from the finite time sin- 
gularity of the original Vieillefosse formulation. The reg- 
ularization is achieved at all scales thanks to the intro- 
duction of an energy-exchange mechanism with smaller 
and larger eddies. A viscous dissipative term is intro- 
duced which acts preferentially at the viscous scale where 
the gradients peak. The model shows very realistic be- 
havior at changing Reynolds numbers, including (i) the 
skewed nature of longitudinal gradients, (ii) the align- 
ment of vorticity with the intermediate eigenvector of 
the strain rate tensor, (iii) the accumulation of events 
along the right tail of the Vieillefosse line and (iv) the 
correct level of intermittency as measured by the Flat- 
ness and its increasing trend at increasing Re\. A set of 



10 




100 1000 Re, 10000 

A, 

FIG. 4: Comparison between a collection of experimental 
data for the longitudinal gradient flatness [2(J, superposed 
with the results from our model (diamonds with error bars). 
Where we have used the variance of the longitudinal coarse 
grained gradient, oa = ((^n 3 ) 2 } 1 '' 2 , to define the Reynolds 
number: Re x = 2E/(3ua A ) 0. Inset: Data for Re x = 1500. 
Normalized Pdf of longitudinal (top) and transverse (bottom) 
band-passed gradients measured at ka- Notice the skewed 
profile for the longitudinal gradients. Longitudinal pdf has 
been shifted along the y-axis for sake of clarity. 

questions remains open. From a dynamical point of view, 
the model seems to under-predict the probability of ob- 
serving strong vortex stretching events - this could point 
to open challenges associated with modeling small-scale 
coherent vortices From a statistical point of view, 
it will be interesting to explore the effects of varying the 
free parameters a and b on intermittency levels. The free 
parameter a defines the relative importance of the local 
Vieillefosse dynamics with respect to the energy exchange 
mechanism. A promising development could be to cou- 
ple the present multi-scale approach with the Lagrangian 
deformation method proposed in (l7j , which proved to be 
successful in reproducing the main features of turbulent 
gradients at small and moderate Reynolds. 
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